Temperature dependent anisotropy of the penetration depth 
and coherence length in MgB 2 
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We report measurements of the temperature dependent anisotropics (7^ and 75) of both the 
London penetration depth A and the upper critical field of MgB2. Data for 7a = A c /A a was 
obtained from measurements of A a and A c on a single crystal sample using a tunnel diode oscillator 
technique. 75 = H^/H^f was deduced from field dependent specific heat measurements on the 
same sample. 7a and 7^ have opposite temperature dependencies, but close to T c tend to a common 
value (7A — 7^ = 1.75 ±0.05). These results are in good agreement with theories accounting for the 
two gap nature of MgB2 . 



The existence of two superconducting gaps in MgB2 |lj 
leads to an unusually strong temperature dependence in 
the anisotropyof the upper critical field, 75 = H^/H^ 
(Refs.|2|,|3j,|J,|5j). The maximum value of 7^ ~ 6 is found 
at the lowest temperatures where it is dominated by the 
anisotropy in the quasi two dimensional a bands. At 
higher temperatures, 7^ is substantially reduced due to 
the increasing contribution of the more isotropic ir band 
5]. The anisotropy of the penetration depth, 7a = A c /A a , 
is predicted to have a markedly different temperature de- 
pendence to that of H C 2 [J E3|- At zero temperature, 
in the clean local limit, the anisotropy in the penetra- 
tion depth is determined only by the anisotropy in the 
Fermi velocity vf and not by the anisotropy in the gap 
• Band structure calculations || show that the average 
vf is approximately isotropic suggesting that A is simi- 
larly isotropic at low temperature^. As temperature is 
increased it has been predicted [J, |fj that 7a increases 
markedly because of the more rapid reduction of the su- 
perfluid density on the it sheets (where the gap is smaller) 
than on the a sheets. 

There have been several experimental measurements 
of the temperature dependence of 7^. Although there 
were some disagreements in early measurements, the 
most recent dataE L 1, agree well with the theoret- 
ical predictions 0, |5| Significantly less experimental 
data has been reported on 7^. Measurements of the dis- 
tortion of the vortex lattice by neutron scattering 
and scanning tunneling spectroscopy suggest that at 
low temperature the anisotropy is small, 7a = 1.2 ± 0.1. 
7a has also been estimated from the anisotropy of H c \ 
ilH cl = / H^f). There has been some disagreement 
between different studies (see Ref. IT2T1 but recent data 
[l2| suggest that at low T, Jh c1 — 1 and there is an 
upward trend in 7jj c1 near T c . However, jH al may not 
be simply related to 7a in MgB 2 because of two gap ef- 
fects. In this paper, we present direct measurements of 
the temperature dependencies of A a and A c in single crys- 



tal samples of MgB2 in the Meissner state, using a sen- 
sitive radio frequency technique. Our data show that 7a 
increases as T — > T c , in agreement with theoretical pre- 
dictions. Measurements of H C 2 anisotropy in the same 
crystals show that the values of 7a and 75 converge as 
T^T C . 

Single crystals of MgB 2 were grown using a high pres- 
sure technique described elsewhere 01 . Some crystals 
used in this study are from the same batch as those 
used in de Haas- van Alphen studies 0], and are hence 
known to be of high quality, with mean free paths of 
4 ~ 850 A, l a ~ 500 A. The T c of the crystals, as deter- 
mined from the midpoint of the heat capacity anomaly, 
was 38.3 K. Measurements of the temperature dependent 
penetration depth were performed using a tunnel diode 
oscillator technique operating at 11.7 MHz jl5j. Samples 
were mounted on a temperature controlled sapphire stage 
(base temperature 1.5 K) with the probe field aligned ei- 
ther into, or along the boron planes. 

Changes in the circuit resonant frequency AF are di- 
rectly proportional to the change in penetration depth 
AA as the sample temperature is varied. With H || c the 
shielding currents flow only in the basal plane, and hence 
the measured frequency shifts are directly proportional 
to AA a . The constant of proportionality is determined 
from the geometry of the sample |lfj| with an accuracy 
of ~ 5 — 10%. For H J_c the shielding currents flow both 
along the c-axis and in the basal plane and the measured 
frequency shifts contain contributions from both AA a and 
AA C . In a rectangular sample |l6| 
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where 2t and 2w are the sample dimensions in the c direc- 
tion and in the in-plane direction (perpendicular to the 
field) respectively. Aivj Lc is the frequency shift observed 
when the sample is withdrawn from the coil, and accounts 
for the demagnetizing factor of the sample. For this study 
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FIG. 1: Temperature dependence of AA a , AA C and AA m i x 
for a single crystal of MgB2. Inset: Ratio of the measured 
frequency shifts for H\\c and HLc. 



samples with simple rectangular shapes and small aspect 
ratios were selected so as to maximize the c-axis contribu- 
tion. The above formula neglects any contribution from 
the top and bottom faces of the sample. For crystals with 
aspect ratios equal to those reported here we expect this 
to be a small correction (< 5%). This is supported by 
the fact that experimentally we find very similar values 
of 7a for samples with different aspect ratios. Most of 
the data presented here are for a crystal with dimensions 
(0.44x0. 33x0. 15)mm 3 with the shortest direction being 
along the c-axis. The accuracy of alignment with the field 
is estimated to be better than 5° in all cases. Corrections 
due to misalignment scale with sin 2 9, and are small [lij . 
The estimated uncertainty in the absolute values of AA a 
is 10%, and 20% for AA C . 

In Fig. ^ we show the measured temperature depen- 
dencies of both AA a and AA C . In this figure we also 
show AA mix = tAF ±c / AF^ C . This quantity is directly 
proportional to the raw frequency shift and would equal 
AA a if the contributions from the c-axis currents were 
negligible. The presence of the AA C contributions en- 
hances AAmix with respect to AA a . We find that the 
temperature dependence of AA C is similar to AA a , but 
a factor ~ 2 larger. This is consistent with the behav- 
ior reported previously 0], although in that study the 
crystals used were not of sufficiently uniform thickness to 
accurately determine A c . 

The penetration depth data are more easily inter- 
preted by calculating the normalized superfluid den- 
sity p — [A(0)/A(T)] 2 . For these calculations we take 
A a (0) = A^(0) = 1000 A in accord with /iSR and neutron 
studies 0, Eil EH ■ We will discuss the effect of varying 
A Q (0) and A c (0) later. 

The calculated superfluid densities p a (T) and p c (T) 
are shown in Fig. Pa(T) and p c (T) have quite 
different temperature dependencies and are both con- 
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FIG. 2: In-plane and interplane superfluid density, p a and p c 
calculated from AA(T), along with fits to the two gap model 
(solid lines). The behavior of p a and p n deduced from the 
fits is also shown (dashed lines). The dotted line shows the 
isotropic BCS weak coupling behavior. Inset: Temperature 
dependent gap functions used in the fit: (dashed and dotted 
lines) taken from Ref. Ilij The BCS gap function (solid line) 
is also shown. 

siderably different to that expected for a conventional 
isotropic BCS superconductor. At low temperature both 
superfluid densities approach unity exponentially, but at 
higher temperature they have a much stronger tempera- 
ture dependence than that expected in the conventional 
case. 

Previously 0], it was shown that p a {T) can be well 
described by a phenomenological two gap model, which 
was first applied to MgP>2 by Bouquet et al. E2 • In this 
model, the contributions to the superfluid density from 
the a and tt bands {p a and p T ) are calculated indepen- 
dently using a temperature dependent energy gap A^ (T) 
which follows the usual weak coupling BCS form but has 
a modified zero temperature value. These two superfluid 
densities are then added to give the total 

p i (T)=Xip*(T) + (l-Xi)p <r (T) (2) 

here i refers to the crystal direction (a or c). The pa- 
rameter Xi sets the relative contribution of the tt band to 
the total normalized superfluid density. Theoretically Xi 
is related to band structure by 

where vf is the i component of vf on the sheet k (k <E 
a, tt) , and the integral is over each pair of Fermi surface 
sheets. In MgB 2 the energy gaps are expected to deviate 
somewhat from the usual BCS T dependence 0, @ . A 
simple correction E3 to the two gap model is to use 
the temperature dependent gap functions as calculated 
using an anisotropic strong coupling model |l8j (see inset 
Fig.HJ. 
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Fits of data to the two-gap model, calculated using 
Afc(T) from Ref. are shown in Fig. [5] As discussed 

previously 0IH, nts to Pa CO S ive A ^ = 29 ± 2 K and 
A CT = 75 ± 5 K, which are in good agreement with values 
obtained from other measurements. The relative weight 
on each band x a = 0.53 ± 0.04, is consistent with Eq. © 
using the calculated Fermi surface parameters 0. 

The strong anisotropy of the a sheets should result 
in the c-axis response being dominated by the more 
isotropic 7r sheets 0. It can be seen in Fig. |21 that p c (T) 
is quite different from p a (T) in a manner consistent with 
the smaller contribution from the a band, and appears 
to be mostly determined by the 7r band where the gap is 
smaller. A fit of the data to the two gap model is shown 
in the figure. Here we have fixed A CT and A„. to the values 
found from the fit to p a (T), and have set x c = 0.91 as this 
produces the correct value for the measured anisotropy, 
7a at T c (see below). The fit does not change appreciable 
if we allow x c to vary by ± 10%, and so is compatible 
which that expected from the band structure (x c =0.99) 
(i.e., the contribution from the a band is not discernable) . 

The fit is slightly worse than for the in-plane data, and 
does not improve markedly if we allow the values of the 
adjustable parameters {x c , A CT and A,r) to vary within 
acceptable limits. We also find the although the exact 
T dependence of A used in the model has little effect on 
the calculations of p a (T) 23], p c (T) and 7a are somewhat 
more sensitive. This suggests that insufficient accuracy 
in the assumed form of A„.(T) are responsible for these 
discrepancies. 

From our measurements of AA and AA C we can cal- 
culate the anisotropy in the penetration depth as a func- 
tion of temperature. This is shown in Fig. IHIa), using 
the same values of A a (0) and A c (0) as used above. Below 
T ~ 7 K, 7a like AA a (T) and AA C (T) is quite temper- 
ature independent, but above this it rises monotonically 
reaching a value of 1.7±0.3 at T c . Changing the value 
of A o (0) from 800 to 1200 A or changing 7A (0) in the 
range 1.0 to 1.2 does not modify the T dependence of 7a 
significantly and also has very little effect on the limiting 
value of 7a at T c [see Fig.|^b)]. The main uncertainty in 
7a (T) comes from the calculation of calibration factors 
used to extract AA a (T) and AA C (T) from the measured 
frequency shifts. To obtain independent estimates of the 
errors we have made measurements on two other crys- 
tals with different aspect ratios (w/t ranges from 8 to 
2.5). These are shown in Fig. Etc). For crystal B, 7a is 
~10% larger at T c but has a very similar T dependence 
to crystal A. Crystal C is again similar but shows more 
upward curvature near T c . We note that close to T c the 
results are particularly sensitive to any macroscopic in- 
homogeneity in the crystal (giving rise to a spread in T c 
values). Data for T/T c > 0.95 are particularly unreliable 
for this reason and are omitted from the figure. 

In Fig. Ota) we show that our data are in good agree- 
ment with the predication of the (parameter free) strong 
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FIG. 3: (a) Penetration depth anisotropy 7a versus temper- 
ature for sample A (symbols). The solid/dashed lines show 
the behavior calculated by Golubov et al. |(| and the phe- 
nomenological two gap model respectively. (b) jx(T) for sam- 
ple A calculated with either 7a(0) = 1.0 (solid symbols) or 
7a(0) = 1.2 (open symbols). The dashed lines show the effect 
of changing A a (0) from 1000 A to 800 A (upper curve) or 
1200 A (lower curve)[with 7a(0) = 1.0]. (c) j x (T) for three 
different samples (A,B,C). 



coupling calculations of Golubov et al. 

la ( 

in the clean 

limit). However, it should be mentioned that Ref. [(| pre- 
dicts A(0) values which are ~ 2 times smaller than the 
experimental ones for samples of the purity we have here. 
We also show in this figure the calculated anisotropy of 
the two gap model based on the fits to the data in Fig. 
121 As mentioned above, we have adjusted x c to give the 
correct value of anisotropy at T c . As for the fits to p c (T), 
the agreement is reasonable. 

For completeness we have also measured the anisotropy 
of H C 2- This was done using both torque and specific 
heat measurements. The heat capacity measurements 
were performed on the same sample (sample A) as that 
used for the A(T) measurements. For the torque mea- 
surements a different, smaller sample was used. The heat 
capacity c p was measured as a function of temperature 
by an a.c. technique [25| in fixed fields up to 7 T ap- 
plied parallel or perpendicular to c. The jump in c p at 
T C (H), was used to deduce H c2 (T) in the two directions. 
In order to calculate 75 (T) , we need to know H C 2 in both 
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FIG. 4: (a) Temperature dependence of 75 and 7a- For 7^ the 
open symbols are specific heat measurements and the closed 
symbols torque measurements. The 7^ data are the same as in 
Fig-EK- (b) Theoreticalpredictions for 7^ (T) and y\(T); solid 
lines - Golubov/Rydh dashed lines - Miranovic/Kogan 

[2^1 . (c) Detail of the experimental data in (a) close to T c . 



field directions at a single temperature. As is almost 
linear with T for H approximately parallel to c we inter- 
polated this data set to give H^iT) at the same temper- 
ature as we have data for H _Lc. Torque was measured 
as a function of H and angle at various temperatures be- 
tween 0.3 K and 35 K using a piezoresistive cantilever 
technique |26| . The angular dependence of H C 2 was then 
fitted to the anisotropic Ginzburg-Landau form to give 
7j . Extracting H C 2 from the torque measurements pill l| 
is more difficult than for heat capacity as pinning effects, 
especially at low temperature, complicate the behavior 
close to H C 2- We find however, that in the region of over- 
lap there is good agreement between the two techniques. 
Our results for 7^(T) shown in Fig. ^ are very similar 
to those reported previously @. In Fig.H we show our 
data for both 75 (T) and 7 a (T). Close to T c it can be 
seen that both anisotropics converge to a single value. 
The c p measurements are particularly well suited to de- 
termining the behavior of 7^ close to T c and we find that 
7 £ (T C ) = 1.75 ±0.05. 

In Fig. we show two different theoretical predic- 
tions for both 7f(T) and 7a (T). The solid lines show the 
results of the strong coupling calculations [7A (T) is as in 
Fig. and 7^ (T) is the theory of Golubov and Koshelev 
|23 usingthe fitting parameters determined by Rydh et 
al. Ref. |9|. The dashed lines show the calculations of 
Miranovic and Kogan |24j (MK) reevaluated for a gap 
ratio Act/Att = 2.7. Both these calculations are in good 
overall agreement with our data. The MK calculation 
underestimates the low T value of 75 because of the as- 
sumed Fermi surface shape (elliptical rather than cylin- 
drical) but correctly predicts that 7^ ~ 7a — 1.9 at T c . 



Although the calculations of Golubov et al. for 7a (T) are 
parameter free, those for 7^ are not and 7^ and 7a are 
only equal at T c for certain parameter values. 

In conclusion, we have measured the temperature de- 
pendent anisotropy of both the penetration depth and 
H C 2 in MgB2. As predicted by theory, they have very 
different temperature dependencies but tend to a com- 
mon value at T c . 

We thank A. Golubov, A. Koshelev, P. Miranovic and 
V. Kogan for giving us details of their calculations. 
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